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I discuss a Bogoliubov inequality for obtaining a rigorous bound on the max- 
imal axial extension of inhomogeneous one- dimensional Bose-Einstein con- 
densates. An explicit upper limit for the aspect ratio of a strongly elongated, 
harmonically trapped Thomas-Fermi condensate is derived. 



The question on the existence and dimension dependence of off-diagonal 
long-range order, for an infinitely extended system of scalar uncharged bo- 
sons, was conclusively answered by Hohenber^. He made use of the fact that 
in the thermodynamic limit of an infinite system of interacting particles of 
bare mass m, the relation (we set ?i = A;b = 1 in what follows, where /cb is 
the Boltzmann constant) 

for the occupation numbers N^. =< > of plane wave states enumer- 

ated by (7^ 0) holds, where T is the temperature and m the bare mass 
of the interacting particles. Angled brackets indicate a thermal ensemble 
(quasi-)averag no and n are the condensate density, associated to = 0, 
and total density, respectively. The inequality is essentially a particular 
representation of the Bogoliubov 1 /k"^ theorem on correlation functionJ^. 

The Hohenberg theorem is based on the mathematical fact that inequal- 
ity leads to a contradiction in dimension D < 2 due to the (infrared) 
divergence of the wave vector integral of which determines the number 
density of non-condensate atoms. We denote the infrared wave vector cutoff 
by kc = 2tt/R, where R is the size of the system, tending to infinity (hence 
kc 0), and the ultraviolet cutoff (the microscopic discretization length) 
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we denote kpi. Taking the integral of both sides of we obtain in the 
one-dimensional case, using N = nR 

The number of excited atoms N — Nq, where A^o is the number of particles 
in the condensate, has to be finite. Therefore, condensates cannot exist at 
any finite temperature, because hq = necessarily in the thermodynamic 
limit of i? ^ oo, due to the divergence of the right-hand side. Physically, 
long-range thermal fluctuations, the influence of which becomes increasingly 
dominant the lower the dimension of the system, destroy the coherence and 
long-range order expressed by the existence of the condensate. 

That there is indeed a difference of trapped Bose-Einstein-condensed va- 
pors of reduced dimensionalit}^^ '^ * ^ * ^ * ^ * ^ * to the thermodynamics in an infinite 
box limit stems from various reasons. First of all, the classification of excited 
non-condensate states by plane waves is not the suitable one in a trapped 
gas: The condensation, in the limit of large total particle number N, takes 
place primarily into a single particle state in co-ordinate spac^^, and con- 
densate and total densities have (in principle arbitrary) spatial dependence, 
n-o no{r), n n{r). Second, there arises the question how interaction 
influences the existence of the condensate, given that inhomogeneity of no(r) 
and n{r). An interacting gas generally has a behavior increasingly different 
from the ideal gas the lower the dimension of the systenP^^ESl, 

For comparison of the long-range coherence properties of trapped Bose- 
Einstein condensates with the spatially homogeneous case of the Hohenberg 
theorem, it is desirable to obtain a relation which is analogous to in as 
close as possible a mannei)^. It is important to stress that a major strength of 
(PQ) is the explicit interaction independence; this holds true as long as this in- 
teraction is velocity independent, i.e., depends only on the (relative) position 
of the particlefl^l. One primary goal thus is to maintain this independence 
of the inequality on particle interactions. 

Like Hohenberg in his original papei^, I start my analysis from the 
Bogoliubov inequality!^, which reads 
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for any two operators A and C. It is valid for any many-body quantum sys- 
tem, for which the thermal (quasi-)averages indicated by the angled brack- 
ets are well-defined and finite. This inequality was used by Mermin and 
Wagner to prove the absence of ferro- and antiferromagnetism in one and 
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two dimensionfl^, in a fashion analogous to Hohenberg's proof of the non- 
existence of off-diagonal long-range order for scalar particles in one and two 
spatial dimensions. 

The operators in relation Q are chosen to be the following smeared 
excitation and total density operators (we keep the spatial dimension D of 
the system for generality): 

A = ifc= / d'^r exp[ik-r]6^{r), C= f d^ r ^l{r) p{r) , (4) 

where T>q is the domain of support of the condensate wave function ^o(r), 
i.e., its (suitably chosen) normalization volume. The excitation operator is 
defined by 6^{r) = $(r) — bo^o{r) and the condensate wave function ^o{r) 
is normalized to unity. This choice of operators is analogous to the original 
Hohenberg theorem derivation irP, where instead of the position space ba- 
sis chosen here, a plane wave basis appropriate for the thermodynamically 
infinite system in a box was taken. In addition, I introduce a variational 
wave number k, which is chosen subject to the constraint that the right- 
hand side of © is maximized (see the example below). In position space, 
the smearing procedure is necessary, because products of two quantum field 
operators diverge if taken at the same point in space. An important step 
in the derivation of a position space version of Hohenberg's theorem is to 
recognize that the so-called mixed response (the "anomalous" commutator) 
is taking the nonlocal form 

p{r), 5d(r')] ) = ^/N^Mr) [-S{r - r') + $*(r)$o(r'')] • (5) 

Here, after carrying out the commutator, < 69 >= ^/Nq =< > has been 
used, where this assignment is valid to 0(l/\/]Vo). The difference to the 
usual Bogoliubov treatment, which would neglect the second term on the 
right-hand side of the canonical commutator [6^{r),S^^ {r')] = 6{r — r') — 
^Q{r)^Q{r') altogether, is that the second term in the mixed response func- 
tion above would disappear in the Bogoliubov approximation. I will argue 
below that neglect of this term in the mixed response, for finite systems, 
would lead to a contradiction with existing, i.e., experimentally realized 
Bose-Einstein condensates. The Bogoliubov approach of [^Oi^q] = vio- 
lates particle number conservation to 0(/), where / = y^l — Nq/N, while 
the present approach resulting in the mixed response ©, violates it only 
to 0(f'^), and this difference, which is quite significant for sufficiently large 
systems, proves indeed crucial for condensates to exisPl. Using the above 
definitions and a position space version of the /-sum rul^^, to eliminate 
details of the Hamiltonian and thus the interaction occurring in the denomi- 
nator on the right-hand side of (jSJ, the final Bogoliubov inequality takes the 
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fornJ^ 

Here, the condensate fraction F = Nq/N and the de Broghe thermal wave- 
length AdB = y2vr/Wr. I define the effective radius of curvature of the con- 
densate to be the curvature radius of the condensate wave function, weighted 
by the total density distribution: 

"1/2 

!^r$o(r-) [-A,$S(r)]n(r)) 
The variational functional C{k) is given by 



no(fc)-$o(fc) / d''r^l{r)\^o{r) 
•JVo 



(8) 



where the Fourier transforms of single particle condensate density and wave 
function are defined to be rio(fc) = / d^r |<l>o('')P 6xp[ifc • r] and ^'o(^) = 
J d^r ^o{r) exp[ik ■ r]. 

There are various notable features of the inequality ©. First of all, it is 
seen that the two-dimensional trapped case is marginal: In the 2D case, the 
radius of curvature Rc scales like Rc oc iV~^/-^. Indeed, in the trapped gas, 
the logarithmic divergence of the 2D case in ((2)) is cut off by the trapping 
potential, and Bose-Einstein condensates can exist even in a (suitably de- 
fined) thermodynamic limil^. Second, the value of C(fc) is strongly reduced 
because of the second term under the square in Q; this is a consequence 
of the fact that we do not make use of the Bogoliubov approximation and, 
therefore, the second term in the mixed response Q is included. Third, 
the relation represents, after the set {Nq, ^Q{r),n{r)} (and thus also F) has 
been specified at a given temperature, an inequality representing primarily 
a geometric statement: It gives a bound on the possible ratio of mean cur- 
vature radius Rc (which is dominated by the weakly confining directions), 
and the de Broglie wavelength AdB- I also mention here that the final term 
in ©, which corresponds to the "—1/2" on the right-hand side of while 
under most circumstances negligible, may be profitably used to obtain upper 
bounds on the possible condensate fraction as a function of temperatur^^. 

To demonstrate that the inequality © is indeed useful, I work out 
an explicit example for an experimentally easy-to-access case. That is, I 
consider a Thomas- Fermi condensate wave function in a quasi- ID harmonic 
trap with lo^ <. which reads ^o{z) = {n^p (l - /A^)^/^. It is 

assumed that F is sufficiently close to unity, A'^g ~ A^, which amounts in 
effect to taking the limit of both sides of © to linear order in 1 — F. The 
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resulting function C is strongly peaked at its global maximum km — 3.7Z^p 
(Am — 1-7Ztf), where C{km) — 1-54 x 10"^. It should be noted that, if one 
were to make the simple "guess" of, say, Am = -^tf, not using a variational 
ansatz for the right-hand side of ©, one would obtain a (much) weaker 
bound, and would therefore strongly limit the usefulness of the inequality. 

Using ((TJ, we have Rc = AZ^p / 3^/2tt , which tell us that, as expected, 
the effective radius of curvature is of order the Thomas-Fermi length of the 
strongly elongated condensate. Neglecting the second term on the right-hand 
side of (inj, at = /cm the rigorous inequality ZTp/AdB < 6 x ^/N — Nq is 
obtained. This result should be compared to the square root of relation 
0, which reads, at — > ra, 

R/^dB < {kpi/n)yJN - No- The bound thus 
obtained on the maximal Thomas-Fermi length of a thin BEC cigar is consis- 
tent with experiments^. It is important to stress, however, that this would 
not be the case, had we neglected, following the Bogoliubov approximation, 
the second term under the square in Q. Doing so reduces the right-hand 
side of the above inequality by about an order of magnitude for quasi-lD 
Thomas-Fermi condensates, and leads to a contradiction with experiment. 
Therefore, in this sense, particle number conservation, violated by the stan- 
dard Bogoliubov prescription, is necessary for the condensate to exist. 

The conditions needed to violate the Thomas-Fermi relation on Ztf be- 
come transparent in experimentally accessible parameters if we write Ztf = 
{6Ndlas/\l - Cas/d±\)^/^{uj±/uj^y/^, where C = 1.4603, a relation valid as 
long as the 3D scattering length is much less than the oscillator length in 
the perpendicular direction, <C dj^. Then 



where / = if all particles are in the condensate, and assumes its maximal 
value / = 1 if there are none. Though the inequality in the above form is 
strictly valid only if / <C 1, for a first estimate of the critical aspect ratio 
it is sufficient to let / = 1. From relation ® it is apparent that, leaving 
iOz (dz), CLs and N fixed, the inequality is most easily violated (apart from 
increasing the temperature), if the aspect ratio is increased by increasing the 
transverse trapping uj±. This route of increasing the aspect ratio to check 
relation (jUj) should in particular be possible in the narrow condensate tubes 
created in 2D optical latticed, where the aspect ratios can be increased to 
1000 and beyond. Indeed, in a such a 2D optical lattice, a transition from 
a superfluid to a Mott insulating phase with no long-range order present 
has been observed, where the point of transition should fulfill inequality 
Also note that close to the geometric scattering resonanc^^ at = 
d±/C, the right-hand side of Q approaches zero and the critical aspect ratio 
strongly decreases, which may be used to detect the scattering resonance. 




d^as 



d± 



(9) 
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To conclude, I stress that the particular advantage of relation © is 
its independence on specific model Hamiltonians for the strongly anisotropi- 
cally confined gas. Bounds for the existence of condensates can thus also be 
given where explicit results are not available, e.g., for long-range interaction 
potentials. By its very nature, as an inequality, it does not (and cannot) 
prove that a certain condensate exists. For such a proof one still has to 
resort to an explicit (many-body) calculation, to show that the state with 
single-particle wave function is macroscopically occupied. However, on 
the other hand, the inequality can rule out parameter ranges for which any 
given set {NQ,^Q{r),n{r)}, at a specified finite temperature, is not possi- 
ble. Essentially, it answers the following concrete question: How much can 
I stretch a strongly elongated, trapped and inhomogeneous Bose-Einstein 
condensate until its long-range coherence breaks? 



REFERENCES 

1. P. C. Hohenberg, Phys. Rev. 158, 383 (1967). 

2. N.N. Bogoliubov, Selected Works, Part II: Quantum and Statistical Mechanics 
(Gordon and Breach, New York, 1991). 

3. A. Gorlitz et al, Phys. Rev. Lett. 87, 130402 (2001). 

4. S. Dettmer et al, Phys. Rev. Lett. 87, 160406 (2001). 

5. H. Ott et al, Phys. Rev. Lett. 87, 230401 (2001). 

6. D. Hellweg et al, Phys. Rev. Lett. 91, 010406 (2003). 

7. H. Moritz et al, Phys. Rev. Lett. 91, 250402 (2003). 

8. T. Stoferle et al, Phys. Rev. Lett. 92, 130403 (2004). 

9. U.R. Fischer, Phys. Rev. Lett. 89, 280402 (2002). 

10. W.J. Mullin, J. Low Temp. Phys. 106, 615 (1997). 

11. E.H. Lieb and W. Liniger, Phys. Rev. 130, 1605 (1963). 

12. D.S. Petrov, M. Holzmann and G.V. Shlyapnikov, Phys. Rev. Lett. 84, 2551 
(2000); D.S. Petrov, G.V. Shlyapnikov and J.T.M. Walraven, Phys. Rev. Lett. 
85, 3745 (2000). 

13. E.H. Lieb, R. Seiringer and J. Yngvason, Phys. Rev. Lett. 91, 150401 (2003). 

14. D. Pines and P. Nozieres, The Theory of Quantum Liquids: Volume I (W. A. 
Benjamin, New York, 1966). 

15. For a particle-number-conserving approach, that is, one which is not breaking 
the global U(l)-symmetry (to a given order in the quantity /), see Y. Castin 
and R. Dum, Phys. Rev. A 57, 3008 (1998). 

16. N.D. Mermin and H. Wagner, Phys. Rev. Lett. 17, 1133 (1966). 

17. E.H. Lieb and R. Seiringer, Phys. Rev. Lett. 88, 170409 (2002). 

18. A.J. Leggett, New J. Phys. 3, 23.1-23.6 (2001). 

19. M. Olshanii, Phys. Rev. Lett. 81, 938 (1998); V. Dunjko, V. Lorent and M. 
Olshanii, Phys. Rev. Lett. 86, 5413 (2001); T. G. Bergeman, M. G. Moore and 
M. Olshanii, Phys. Rev. Lett. 91, 163201 (2003). 



